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' The instabilities induced on a two-dimensional system of correlated electrons by the anisotropies 

pH ^ of its Fermi line are analyzed on general grounds. Simple scaling arguments allow to predict the 

d • opening of a superconducting gap with a well-defined symmetry prescribed by the geometry of the 

Fermi line. The same arguments predict a critical dimension of 3/2 for the transition of the two- 
J — , dimensional system to non-Fermi liquid behavior. The methods are applied to the t — t' Hubbard 

^SJ ' model in a wide range of dopings. 
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I. INTRODUCTION 



Two-dimensional electron systems are the stage in which most relevant quantum phase transitions arise. The 
prototype of these phenomena is the Mott-Hubbard transition , which can only be understood in the context of 
' ^ ' a very strong correlation between the electrons. Many of the properties of the high-Tj, superconductors are also 
Q , supposed to be the consequence of strong correlations in the two-dimensional layers In both cases, these reflect 
Q themselves in the presence of very marked features in the single-particle dispersion relation. The case of the copper- 
. oxide compounds has been actually studied with increasing accuracy from the experimental point of view, and the 
I ■ angle-resolved photoemission measurements have been able to characterize the drastic change that the dispersion 
^ , relation suffers when the stechiometric material is increasingly doped In the antiferromagnetic insulator regime, 
^\ • the characteristic peaks at {tt/2, 7t/2) and symmetry related points are a consequence of the doubling of the unit cell 
0^ , and the hybridization of modes by the wavevector (7r,7r) [Q. In the superconducting regime, the two-dimensional 

■ layers tend to develop very flat bands at the points (7r,0) and (0, tt) |||,|^. It is not well-known the mechanism by 
, which the electron system undergoes such an strong renormalization but, reversing the line of reasoning, one may ask 

■ to what extent the pronounced features of the Fermi surface are at the origin of the unconventional properties in the 
, superconducting regime. More generally, there is a lack of some dynamical formulation predicting the evolution of the 
' Fermi surface upon doping and, at this point, it is only possible to characterize the shapes that lead to the opening 
I of a more stable phase (Mott insulator, superconductor). 

^ In general, phases different from the normal metallic state may arise from instabilities that have their origin in 

^ the peculiar geometry of the Fermi surface. Such instabilities may lead in some instances to antiferromagnetism, 
(— ( superconductivity, or the formation of a charge-density-wave structure. A well-known example is the case of a Fermi 
Q ] surface with the property of nesting, that is, with two portions that are mapped one into the other by some fixed 
O momentum Q. The system shows then an enhanced response to a perturbation with such wavevector and, if the 
^ ' corresponding modulation is commensurate with the lattice, it may lead to antiferromagnetic order or to a charge- 
density- wave structure, depending on the character of the interaction Q. On the other hand, there are also instances 
. , in which the geometry of the Fermi surface may give rise to a superconducting instability, starting from a pure 
' repulsive interaction. The basis of this mechanism was laid down by Kohn and Luttinger back in 1965 In the 

case of a 3D electron system with isotropic Fermi surface, there is an enhanced scattering at momentum transfer 
2kF, which translates into a modulation of the effective interaction potential V{r) ^ cos{2kpr)/r^. This oscillating 
behavior makes possible the existence of attractive channels, labelled by the angular momentum quantum number. 
The superconducting instability arises from the channel with the strongest attractive coupling. In two dimensions, 
however, it has been shown for the weakly non-ideal Fermi gas that the effective interaction vertex computed to second 
order in perturbation theory is independent of the momentum of particles at the Fermi surface | [lO|] . In this case one 
has to invoke higher-order effects in the scattering amplitude in order to find any attractive channel [ pT| . Otherwise, 
superconductivity due to the Kohn-Luttinger mechanism also arises in 2D systems with a Fermi surface that deviates 
from perfect isotropy. It has been studied in the Hubbard model at very low fillings, for instance, with the result 
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that the dominant instabihty may correspond to d^y, d^2_y2 or p-wave symmetry, depending on the next-to- nearest 
neighbor hopping Jl^-p^. 

In the present paper we review the possibihty of Kohn-Luttinger superconductivity in the case of Fermi surfaces 
with a shape appropriate to the description of the two-dimensional hole-doped copper oxide layers. As we want to 
consider the regime around the optimum doping for superconductivity, this implies to deal with highly anisotropic 
Fermi surfaces. We recall that a most accurate description of the electronic interaction in the layers of the hole-doped 
cuprates is given by the Hubbard model with one-site repulsive interaction and significant next-to-nearest-neighbor 
hopping i' < jl^ . The values that give the best fit to the Fermi lines determined by the photoemission experiments 
seem to be around —0.3 times the t hopping parameter. Typical energy contour lines for that value are shown in 
Fig. I. There we can appreciate how the topology of the Fermi line changes upon doping after passing through the 
saddle points at (7r,0) and (0, tt). Quite remarkably, for filling levels right below the Van Hove singularity the Fermi 
line develops inflection points, i.e. points at which the curvature changes sign. The presence of these points leads 
to a strong modulation of the effective interaction of particles at the Fermi line, as their scattering is enhanced at 
momentum transfer connecting every two opposite inflection points. This effect is reminiscent of what happens in the 
case of nesting of the Fermi surface, although in the present circumstance there are no finite portions of the Fermi 
line that are mapped by a given wavevector. As a consequence of that, there is not a marked tendency towards a 
magnetic instability, but the dominant instability corresponds to superconductivity by effect of the modulation along 
the Fermi line of the effective interaction between particles with opposite momentum fl^ . 



(0,7t) 




FIG. 1. Different shapes of the Fermi line for the t ~ t' Hubbard model about the Van Hove fiUing. Two opposite inflection 
points are marked on the figure. 



We will see that, in the mentioned cases of strong anisotropy of the Fermi line, there is a direct correspondence 
between its topology and the channel in which the dominant superconducting instability opens up. As we will be 
carrying out the discussion for the t — t' Hubbard model, the channels are characterized by the irreducible represen- 
tations of the discrete symmetry group. We will show that, in the case in which the Fermi line has a set of evenly 
distributed inflection points, superconductivity due to the Kohn-Luttinger mechanism takes place in the extendend 
s channel, with nodes at the inflection points. On the other hand, when the Fermi line is above the saddle-points, 
the superconducting instability most likely takes place in the dx2^y2 channel, with nodes at the points where the 
curvature reaches a minimum. Between these two instances, we find the situation in which the Fermi energy sits at 
the Van Hove singularity, what requires special care in order to deal with the divergences that appear in perturba- 
tion theory associated to the logarithmic density of states. In the context of highly anisotropic Fermi surfaces, the 
implementation of the Kohn-Luttinger mechanism was proposed within the investigation of the electron system near 
Van Hove singularities (it). Related analyses have been also undertaken by other authors []l8|-^. As we will see, 
such an electronic mechanism of superconductivity near a Van Hove singularity may be relevant to the physics of 
the high- Tc cuprates [2l|-|2^]. From the technical point of view, this is due to the fact that the modulation of the 
effective interaction vertex is driven then by a marginal operator, rather than by an irrelevant operator as it happens 
in general, so that in some region of the phase diagram critical temperatures may be reached as high as those found 
in the cuprates [ p4| . 

We will adopt a renormalization group (RG) approach to the discussion of the subject. In recent years, RG methods 
have been applied to the description of interacting electron systems [0,p5[ . They have given a precise characterization 
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of Fermi liquid theory, as well as an efficient classification of its relevant perturbations Q . Among a reduced number 
of them, superconductivity appears as one of the possible instabilities. It only requires the presence of an attractive 
coupling in some channel, though small it may be in the bare theory. RG methods are also particularly well-suited 
to the discussion of highly anisotropic Fermi surfaces, since a change of the topology is accompanied by a change of 
the scaling dimension of the interaction vertex. This is the way in which the enhancement of scattering produced 
by some features, as for instance the inflection points, is understood in the RG framework The RG approach 
is specially useful when dealing with unconventional mechanisms of superconductivity, as it may face the possibility 
that the normal state does not fall into the Fermi liquid description. Under this circumstance, it may not be justified 
to use a perturbative approach, or even to rely on certain sum of diagrams in perturbation theory, to study the 
superconducting instability. Yet in the RG approach the only basic assumption is that there is a well-defined scaling 
near the Fermi level, which allows to identify the scaling operators in the low-energy effective theory. Apart from 
addressing the question of superconductivity, the discussion can be devoted to establish when such scaling conveys 
to non- Fermi liquid behavior , spoiling the conventional picture that considers superconductivity as a low-energy 
instability of Fermi liquid theory. 

The paper is organized as follows. In the next section we briefly review the renormalization group approach to 
interacting electron systems. In section 3 we investigate under which conditions non-Fcrmi liquid behavior may arise 
by effect of the anisotropy of the Fermi surface. The analysis of the superconducting instabilities in the presence of 
inflection points in the Fermi line is carried out in section 4. The RG approach is adapted to the discussion of the 2D 
electron system near Van Hove singularities in section 5. Finally, our conclusions and outlook are drawn in section 6. 



II. RENORMALIZATION GROUP APPROACH TO INTERACTING ELECTRONS 

The wilsonian RG approach, that has proven to be so useful in the study of classical statistical systems, has been 
recently implemented in the investigation of many-body systems with a Fermi surface [^. In this context, the basic 
idea of the method is closely related to the concept of effective field theory |^ . One is interested in identifying the 
elementary fields and excitations at a very low energy scale about the Fermi level, that is, at a much smaller scale 
than the typical electron energies of the bare theory. The strategy to accomplish this task is to perform a progressive 
integration of high-energy modes living in two thin shells at distance A in energy below and above the Fermi surface. 
Sufficiently close to it, the fields of the effective theory have to scale appropriately under a reduction of the cutoff 
A —f sA. Then, one can make an inspection of all the possible terms built out of them contributing to the effective 
action of the theory, checking their behavior under the scale transformation 

gi [ dtd'^p ^ gi s"' I dtd^p Oi (1) 



The terms in the effective action that scale with an exponent > are said to be irrelevant, as their effect becomes 
weaker and weaker close to the Fermi surface. On the other hand, if there appear some terms scaling with ai < 0, 
this means that we have chosen a wrong starting point for the effective field theory, since there are couplings that are 
not stabilized already at the classical level. The interesting case corresponds to having all the exponents > and, 
in particular, operators with = (so called marginal operators), as this is the situation where the RG approach 
can address the existence of a fixed-point of the scale transformations in the low-energy theory. 

The system of interacting electrons with a sphere-like Fermi surface provides a good example of how the above 
program is at work [ p7[ . We focus from now on in two spatial dimensions and write the action with the most general 
four-fermion interaction term 

S = j dtd^p (z*+(p)at*,(p)-(e(p)-£F) 

+ j didVrf^P2d^P3rfVC^(Pi,P2,P3,P4) 'J'+(Pi)1'J-(P2)*^'(P4)*<t(P3)'^(P1+P2-P3-P4) (2) 

We suppose that the integration of high-energy modes has already proceed to energies close to the Fermi level, so that 
^'(p) stands here for some renormalized scaling fermion field. The important point to notice is that, at each point 
of the Fermi line, only the component of the momentum orthogonal to it scales with the energy cutoff. In fact, near 
the Fermi level we can decompose every momentum p into a vector P to the closest point in the Fermi line and the 
orthogonal component p^, as shown in Fig. |2[ 

P = P + P± (3) 
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Upon the scale reduction of the energy cutoff, we have p± — > sp±. 




FIG. 2. Decomposition of the momentum p into a transverse and a large component at the Fermi line. 

It becomes clear, for instance, that in the noninteracting theory the field ^'(p) has a well-defined scaling rule that 
makes the first term of scale invariant. We can rewrite the free effective action in the form 

^0 ~ y dtdPdp^ {i^t{p)dt^Ap) - «f(PW*+(p)*.(p)) (4) 

Under a change of the cutoff A —^ sA and corresponding changes in the scaling variables t — > s^^t,p± — > sp±, the 
field has to transform according to ^(p) — > s^^/^5'(p) in order to keep (||) marginal. 
A similar analysis applied to the interaction term of the effective action 

Sint ~ J dtdPidpij_dP2dp2±dP3dp3±dPidp4j_ 

C^(Pl,P2,P3,P4) «'^(Pl)*J'(P2)1'<T'(P4)«'^(P3)'5(pi +P2 -P3 -P4) (5) 

leads to the conclusion that, in general, it scales in the form Smt — *■ sSmt- We are assuming that, upon scaling, 
the orthogonal component of the momenta are irrelevant in the momentum-conservation delta function, compared to 
the large P components, so that 6{pi -t- P2 — Pa — P4) ~ S{Pi -I- P2 — P3 — P4). Thus the interaction turns out to 
be irrelevant for generic values of the momenta (assuming implicitly a smooth dependence on them of the potential 

C^(Pl,P2,P3,P4))- 

The important remark put forward in Refs. |^ and psj is that there are special processes in which the kinematics 
forces a different scaling behavior through the constraint of momentum conservation. This is the case when the 
combination of the momenta at the Fermi surface identically vanishes. Pi -I- P2 — P3 — P4 = 0. When this happens, 
the scaling of the orthogonal components of the momenta is transferred to the scaling of the delta function, and 
the four-fermion interaction term becomes marginal. This explains why in Fermi liquid theory there are only a few 
channels that are not irrelevant in the low-energy theory. 

Specifically, the above identity for the sum of the components at the Fermi surface is satisfied by i) Pi = — P2 
and P3 = — P4, what characterizes the so called BCS channel, ii) Pi = P3 and P2 = P4, that leads to the forward 
scattering channel, and iii) Pi = P4 and P2 = P3, that differs from the previous one by the exchange of the outgoing 
particles. In this classification we have not taken into account the spin of the particles, but it is clear that, if in the 
third case the spin of P4 is different to that of Pi, the scattering process is fully differenciated from that of forward 
scattering. In the case of a spin-dependent interaction, it makes sense to consider iii) as a different channel on its 
own, that we will call exchange scattering channel |2^. The three different channels leading to marginal interactions 
are represented graphically in Fig. ^. 
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FIG. 3. Marginal interactions in Fermi liquid theory. 

We address now the behavior of the marginal interactions under quantum corrections by focusing again on a system 
with a Fermi line having the topology of the circle. This is actually the case best studied in the literature. We borrow 
the main conclusions from Refs. Q and p5[ |, which state, to begin with, that the interaction in the forward scattering 
channel remains unrenormalized in the quantum theory. That is, the integration of the high-energy modes does not 
produce any dependence of such marginal interaction on the cutoff by effect of the loop corrections. This is the way 
in which the Landau theory of the Fermi liquid is recovered in the present context, with a finite correction of the bare 
parameters and a F function that encodes all the information about the four-fermion interaction. 

On the other hand, the pairing interaction in the BCS channel gets dependence on the cutoff by integration of 
modes in the two thin slices at energies A and —A. At the one-loop level, for instance, the particle-particle bubble is 
the diagram that contributes to the renormalization since, if one of the particles in the loop carries momentum p at 
the slice being integrated, the other has a momentum — p that also belongs to the slice. In two dimensions, taking 
into account the dependence of the BCS coupling on the angles 9 and 9' of the respective incoming and outgoing 
pairs, the one-loop differential correction becomes 

VA+dAi9,9')^VAi9,9') + N{eF)dA j d9"VK{9,9")VK{9" ,9') , dA < (6) 

where N{£f) is the density of states at the Fermi level. 
The one-loop flow equation 

p^^X^Ml = Niep) J d9"V{9, 9")V{9", 9') (7) 

describes the approach to the Fermi line as A ^ 0. Thus, in systems where the original bare interaction is repulsive 
the BCS channel has a general tendency to become suppressed at low energies. It is only when V or any of its normal 
modes is negative that the renormalized coupling becomes increasingly attractive near the Fermi level. This is the way 
in which the superconducting instability arises in the RG framework [^^. The main conclusion is that the system 
with uniform repulsive interaction and isotropic Fermi line provides a paradigm of the Landau theory of the Fermi 
liquid, since all the interactions turn out to be irrelevant but the F interaction, that remains finite and non-zero close 
to the Fermi level. 

The above RG scheme gives also the clue of why the slight deviation from isotropy may trigger the onset of 
superconductivity. In general the analysis has to be carried over the couplings Vn for the normal modes appropriate 
to the symmetry of the Fermi line. The flow equation may be decomposed then into 



= N{eF)V^ (8) 



and the flow of the couplings is given by 



' l + N{eF)Vn{K) log(Ao/A) ^ ' 

If, at some intermediate scale Aq, a coupling V^i(Ao) becomes negative, then a superconducting instability develops 
in the corresponding mode. The Kohn-Luttinger mechanism is based on the fact that this happens in practice, 
even for an isotropic Fermi surface in three dimensions, since the corrections to the bare V{9, 9') by particle-hole 
diagrams already introduce the degree of anisotropy needed to turn some of the Vn negative [Q. In two dimensions, 
the anisotropy of the Fermi line is also a necessary condition to produce an attractive interaction in some of the 
normal modes, as pointed out in Ref. ||l^, imless the particle-hole corrections to the bare V{9,9') are carried out 
beyond second order in perturbation theory pT| . 
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Previous analyses of the Kohn-Luttinger mechanism have focused on small deviations from isotropy of the Fermi line 
1^ . The relevance of the superconducting transition is related anyhow to the source of anisotropy that is producing 
the attractive coupling in the model. An attractive but very tiny T^„(Ao) may lead to superconductivity at a very 
small scale E ~ Aoexp{l/ (iV(eF)F„(Ao))}. Our RG analysis makes clear that the anisotropy has an effect through 
corrections to the bare scattering V{6,6') that are irrelevant as the cutoff A 0. We will see that they scale in 
general like A^/^, that is, with the same scaling law as for scattering with 2fci? momentum transfer in the isotropic 
Fermi line. On the other hand, we will be dealing with topological features of the Fermi line, i.e. inflection points, 
that enhance the attractive modes to scale like A^/"* Inflection points can be thought of as precursors of the 

development of the saddle points at the Fermi line, a situation in which the particle-hole corrections to scattering 
processes become marginal A" ). Consequently, in the context we are facing of a repulsive bare interaction, the 
superconducting instability is significantly strengthened near the Van Hove singularity p"7j~ p0| , p9[ , although then there 
is also competition with magnetic instabilities in the system p^ , pO| - p2| , p4| . The RG analysis of the problem becomes 
quite delicate, moreover, as one has to derive a proper scaling in a model which has enhanced log^(A) divergences in 
the BCS channel. The special topology of the Fermi line passing by the saddle points introduces significant changes 
in the analysis carried out for the Fermi liquid theory and it requires the specific treatment that we will discuss in 
Section 5. 

III. INFLECTION POINTS OF THE FERMI LINE 

As discussed in the preceding section, the scaling of the different couplings can be read off from the cutoff dependence 
of the particle-particle and particle-hole diagrams which modify them. When the cutoff is sufficiently close to the 
Fermi surface, these diagrams are reduced to convolutions of single particle propagators taken at the Fermi surface. 
The relevance or irrelavance of the couplings is determined by the density of states of electron hole pairs, or of two 
electron states. 

We can deduce the real part of a given diagram from its imaginary part by performing a Hilbert transform. The 
imaginary parts measure directly the density of excited states, and are usually easier to estimate. Let us now imagine 
that we can determine that, for a given susceptibility, x, the low energy behavior is Imx(q, w) ~ cj". Then, Rex(q, a;) 
includes a term with the same lu dependence if a 0, and shows a logarithmic behavior if a = 0. This contribution 
will dominate the corrections to the couplings when it is present. 

The quantity Imx(q, w) is a measure of the electron-hole or electron-electron excitations of total momentum q and 
energy u. The simplest way to estimate it is to measure the total number of excitations with momentum q and energy 
equal or less than lu. By deriving the volume of phase space which fulfills these contraints with respect to w, one 
obtains Imx(q, w). 

The leading diagram which has a non trivial low energy behavior is the electron electron propagator at zero 
momentum, the BCS channel, as discussed in the previous section. It is easy to show that ImxBCsil ^ 0,llj) ^ 
constant as a; — *■ 0. This diagram describes the main instability of Fermi liquids. In the following, we analyze the 
modifications that insertions of electron-hole propagators induce in this channel, following the ideas of Kohn and 
Luttinger. 

We first need to estimate the the energy dependence of the electron hole propagator Imx(q, a;). In one dimension, 
it is relatively straightforward to show that the imaginary part of the electron-hole susceptibility, Imx((7, w), is finite 
all the way to zero energy only when q — 2kp. The limit of lmx{2kF,oj) as goes as the single particle density 
of states at the Fermi level (a schematic view of the available density of states is shown in Fig. ^ ). 

From limi^_,o Inix(2fc_F, oc ImG(ei?), we deduce that linii^^o Rex(2fci?, w) cx ImG(ei?) log(A/aj), where A is the 
cutoff of the theory. Outside q = 2kp there are no gapless electron- hole excitations in one dimension. The log 
dependence in Rex, when inserted in the appropiate diagrams, leads to the existence of marginal couplings, and to 
the non trivial phenomenology of Luttinger liquids. 

In higher dimensions, Imx(q, tj) remains finite at low energies when q connects two points at the Fermi surface. 
The density of electron-hole pairs at this wavevector depends on the curvature of the Fermi surface at these two 
points, and also on their relative orientation. In general, the two patches of the Fermi surface are not parallel. We 
can linearize the dispersion relation as e(ki -|- k') = v^k', where ki — ko and k2 = ko + q denote two points at 
the Fermi surface. The constraints described earlier imply that: w > e(ki -I- k') — £(k2 + k'), and, simultaneously, 
e(ki -I- k') > and £(k2 -I- k') < 0. These three conditions, if -v^. and v|, are not parallel, define a triangle in k' space 
such that the basis and the height are proportional to w, as schematically shown in Fig. |^. Thus, Imx(q, t^^) oc lo. 
This linear dependence on lo remains when the local density of electron-hole excitations is determined, by integrating 
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over q. This is the relevant quantity which determines the couphng of impurities to the electron liquid, and can lead 
to non trivial divergences in perturbation theory, as best shown in the Kondo model. 

The previous analysis does not hold when q connects two regions of the Brillouin Zone where the Fermi surface is 
parallel. The simplest case is also shown in Fig. |[ The linear expansion of the dispersion relation is not enough to 
estimate Im^. Using a frame of reference such that one axis is perpendicular to the two Fermi surfaces, one can write: 

e(ki +k') = v],k'^+aikf 

e(k2 + k') = vj-k'^ + aafcf (10) 

The boundaries in k' of the region where particle-hole excitations of energy less than oj are possible are given by 
k'^ ~ io/{vp + Vp) and k'n ~ y/uj/{ai + 02), giving a total area oc u)^^"^. Hence, Imx(q, t^^) ~ v^- This analysis 
is restricted to two spatial dimensions. In higher dimensions, one has to take into account the additional D ~ 2 
dimensions. Along each of them, the available area is limited to a vector in k space which also scales as ^/lo. Hence, 
the volume is oc and Imx(q, cj) '-^ There is a logarithmic divergence in Rex when D = \, which 

plays the role of the upper critical dimension. 

We now assume that the Fermi surface is sufficiently anisotropic, as shown in Fig. ^, so that there are inflections 
points at wavevectors {k^}. A momentum transfer equal to q^ = 2ki connects the two inflection points at ±ki which, 
in addition, define parallel patches of the Fermi surface. The expansion in Eq. ([l0| ) is not enough, as ai — a2 — 0. It 
must be replaced by: 

e(ki + k') = WFfc± + /3fcf + 7fcf 

e(k2 + k') = -vpk'^ - Pkf + 7A:f (11) 

where we are using the symmetry between the two points. The area with electron-hole pairs with energy equal or less 
than w is now limited by k'j_ ^ uj/vp and k'^^ ~ {uj/^Y/^. Hence, Imx(q, w) ~ w^/''. 

The previous argument can be extended to arbitrary dimensions. In general, an inflection point along one of the 
D — 1 transverse directions at the Fermi surface will not coincide with an inflection point along another direction. 
Thus, in the remaining D — 2 directions, the phase space is limited by quadratic terms, obtained from an expansion 
like that in Eq. ([lO|). Hence, Imx(q, w) ^ tj(^-2)/2+i/4_ rpj^^ dependence of Rex on the high energy cutoff of the 
theory becomes logarithmic when D — 'i/2 ||l^]. In the RG language, the couplings at wavevector q are marginal, and 
require a non trivial renormalization, at this dimension. It is interesting to note that, for isotropic Fermi surfaces, the 
critical dimension is 1 ]3^j3^ . 



IV. ANISOTROPIC FERMI LINES AND KOHN-LUTTINGER SUPERCONDUCTIVITY 

In the previous Section we have seen that the strength of electron scattering bears a direct relation with the 
curvature at each point of the Fermi line. When this becomes highly anisotropic, the strong modulation in the 
angular dependence of the scattering between particles with momentum p and — p may lead to a superconducting 
instability, as explained in Section 2. We will address here the effect of inflection points in the t — t' Hubbard model, 
that is the simplest model with this kind of features in the Fermi line, and we will determine the symmetry channels 
in which pairing can take place. 

The main issue is whether there exists any attractive coupling, at some intermediate energy scale Aq, for some of 
the normal modes of the bare BCS vertex V{Kq). Here Aq is supposed to be a small fraction of the whole bandwidth, 
but not yet at the flnal stage of the renormalization process. We recall that V depends on the angles 9 and 9' of 
the respective incoming and outgoing pairs of electrons, so that it may be decomposed into eigenfunctions with well- 
defined transformation properties under the action of the lattice symmetry group D^. This has four one-dimensional 
representations, labelled respectively by Ai,A2, Bi and B2, and a two-dimensional representation labelled by E. The 
complete sets of eigenfunctions for the respective irreducible representations are given by 
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Thus, the BCS vertex may be expanded in the form 

V{e, e') ^Vq + Vi {cos{e) coa{e') + sm{e) sm{9')) + V2 cos(26l) 008(26*') + V3 am{2e) sm{2e') 

+Vi (cos(36i) cos(36'') + sm(36') sin(36'')) + V5 coa{4e) coa{4e') + Vq sin(46') sin(46'') + . . . (13) 

According to the discussion in Section 2, it suffices that any of the Vn couphngs becomes negative for a superconducting 
instabihty to take place in the process of renormaUzation, with the corresponding symmetry of the order parameter 
— s-wave in the case of Ai and A2, and d-wave in the case of the Bi and B2 representations. 

We know that aU the interactions are irrelevant in Fermi liquid theory, except the marginal V, F and E interactions. 
Furthermore, a repulsive interaction in the V channel flows to zero upon renormalization towards the Fermi level. 
The point is that, even starting with a bare repulsive interaction as in the t — t' Hubbard model, some of the irrelevant 
couplings may lead to attraction in the BCS channel well before arriving at the final stage of the renormalization 
process. Under these circumstances, one needs to make an inspection of the theory at the intermediate scale Aq, in 
which the irrelevant interactions discussed in Section 2 have not become negligible yet. They may turn negative, in 
fact, some of the Vn couplings. When this happens, by considering Aq as the starting point of the renormalization 
group flow one observes the appearance of a strong attraction that leads to superconductivity in the low-energy 
effective theory. 

Before switching on the quantum corrections, the marginal interaction V equals the bare on-site repulsion U of the 
t — t' Hubbard model. Given that the corrections we want to study at the intermediate scale Aq are given by irrelevant 
operators, we may estimate their effect in perturbation theory. At the one-loop order, the diagrams that contribute 
in general are those in Fig. |^. In the case of the local Hubbard interaction, we take the convention of writing it as 
mediated by a potential between currents of opposite spin. This makes clear that the first two diagrams in Fig. |^ 
cannot be drawn with the interaction at hand, since they require an interaction between currents with parallel spins. 
They are present in the case of more general short-range interactions, but they tend to cancel out when such range 
shrinks to zero. 

The contribution Il{9,9') of diagram (c) has the overall effect of reinforcing the bare repulsion between the elec- 
trons. We recall that and 0' label the angles between the incoming momentum p and the outgoing momentum k, 
respectively, and the x axis. When the two angles are the same, the momentum flowing in the interactions is large 
and connects opposite points at the Fermi line. We may borrow then the results of the previous Section to determine 
the scaling of the particle-hole contribution ll{0,0) correcting V{0,d). 

The simplest case to discuss is that of a set of eight inflection points on the Fermi line evenly distributed in the 
angular variable 0. Then the function ll{0, 0) must admit an approximate representation of the form 

\l{0,0)f=i a + bs\T?{'i0) (14) 

We know from the previous Section that the coefficients a and b have different order of magnitude, given their different 
dependence on the energy cutoff Aq. At = 0, the scattering does not differ much from that of momentum transfer 
2kF on an isotropic Fermi line, so that a ~ 0{h^^'^ / 0^^'^)^ with the notation of the previous Section. At 6* = 7r/8, the 
scattering is greatly enhanced since it takes place between particles sitting at opposite inflection points on the Fermi 
line, implying 6 ~ 0{h}J^ J^^^'^). 

The approximate expression for the bare BCS vertex at the energy scale Aq 

U + TI{0,0)kU + a + hsm^{A0) (15) 

matches well the expansion (^3|) truncated to include up to the Ve term. This allows to find the scaling of some 
combinations of the coefficients 

Vo + Vi + V2 + Vi + V^-U ^ 0[Kl'^) 
V:i^V2^Q 

V^-V^^ 0{Kl") (16) 

In order to close this system of equations we need an additional piece of information, that we get by estimating 
Ii{0,—0). This object stands for the scattering with a momentum transfer connecting a point on the Fermi line 
and its symmetric with respect to the y axis. Therefore, it has again the characteristic dependence ~ ©(Ag^^/a^/^) 
if 6' = or TT and it shows a crossover to a behavior ~ 0(Ao) for other values of 0. We may approximate then 
n(6', —0) Ki c + dco&^{0), where c ~ 0(Ao) and d ~ 0(Aq''^). Comparing with (|l^), this imphes that 
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Vo-Vx-V:i-Vi-V^~U ^ 0(Ao) 
2^1 ^ O(Ay') 

1/2 + ^3 ^ 
^4 

^5 + - (17) 

Putting together (|l6|) and (jl^), we obtain Vq ^ -V^ ^ Oikl''^), Vq-U ^ Oikl''^) and Vi - O(Ay^), the rest of the 
coefficients being higher-order powers of Aq^^. 

Our analysis shows that there is a negative coupHng at an intermediate energy scale Aq in the expansion in 
normal modes of the marginal interaction V{9,0'). The negative coupling V5 sets actually the leading behavior of 
the coefficients Vm for n ^ 0, as the energy cutoff is sent to the Fermi level. The approximations we have made of 
Ii{8,6) and 11(6', —0) by specific periodic functions do not have therefore major influence in the determination of the 
attractive channel, as long as the higher harmonics in the expression (|l^) can be considered subdominant with respect 
to the behavior of the couplings Vi, Vi; and Vq. 

We conclude that there must be a range of dopings, around the filling level where the inflection points are evenly 
distributed in the angular variable on the Fermi line, in which a superconducting instability opens up in the channel 
corresponding to the Ai representation of the lattice symmetry group |p^ . This means that the order parameter 
of superconductivity has the so-called extended s-wave symmetry, with nodes at the inflection points of the Fermi 
line. This location of the nodes has to be shared to a certain degree of approximation by any Fermi line with a set 
of inflection points. The gap in the superconductor tends to close up at the points where the interaction gets more 
repulsive. In the model with local interaction, this happens at the points where there is more phase space for the 
scattering between the particles. The character of the bare interaction is crucial in this consideration, since it relies on 
the fact that the main effect of the anisotropy at the scale Aq comes through diagram (c) in Fig. This contribution 
goes in the opposite direction to screening the bare interaction. It becomes clear that, under a more general type of 
interaction, different possibilities for the opening of an attractive channel could arise, while the methods outlined in 
this Section should still be useful to study very anisotropic Fermi lines. 

In the context of the t ~ t' Hubbard model, the situation in which the Fermi line approaches the saddle points 
at (7r,0) and (0,7r) has also great phenomenological interest. One possibility would be to address this problem by 
means of a sequence of different dopings and Fermi lines having the kind of inflection points we have just discussed. 
However, this would imply that, when approaching the saddle points, two inflection points would become quite close 
near (tt, 0) and (0, tt). It is clear that the kind of approximations we have made before should break down at a certain 
point in this limit, as the distribution of the strength of the scattering along the Fermi line becomes quite sharp. In 
the limit case where the Fermi level is at the Van Hove singularity, the distribution of the curvature of the Fermi line 
becomes actually singular at the saddle points. This suggests that a different method has to be deviced to deal with 
this special situation, on which we will elaborate in the following Section. 

Anyhow we still have the alternative of approaching the saddle points from above the Van Hove singularity, with 
the kind of round shaped Fermi lines shown in Fig. |l|. As far as the curvature is distributed smoothly over the Fermi 
line, we can approximate again the ECS vertex by a certain number of harmonics, the rest of them being much more 
irrelevant as the cutoff is reduced. 

For the round Fermi lines depicted in Fig. Qthe bare vertex H(0, 9) is ~ 0(Ay^/a^/^), so that it is just modulated 
by the local curvature a of the Fermi line. It is a function with period equal to 7r/2, reaching maxima at the angles 
where the curvature has the minimum value. For convenience, we measure now angles with the origin at (7r,7r). We 
may take then 

Tl{e,0)7ii a + hs\r? {29) (18) 
where a ~ 6 ^ 0{h^^^). Comparing with the expansion (^3|), we obtain the relations 

V0+V1+V2-U ^ o{Ay^) 

V^-V2^ 0{Ky^) (19) 

Using as before the estimate T[{9,—9) « c -I- dcos^{9), with c ^ 0{Kq) and d ^ ©(A^^), we have the additional 
constraints 
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Vo-Vi-Vs-U ^ 0(Ao) 

V2 + V3^0 (20) 

Taking into account ^ and (|o|), we conclude that V3 ¥2-^ 0{AI^^) and Vi - ©(A^^). In this case, we find 

that among the dominant contributions as Aq — > there is a negative couphng, V2, that corresponds to an instabihty 
in the d^^-y^ symmetry channel. This result concerning the order parameter is in agreement with more detailed 



calculations for Fermi lines of similar shapes 1 20 



Thus, although we are not able to deal yet with the particular instance in which the Fermi level sits at the Van Hove 
singularity, the above discussion makes plausible that a strong superconducting instability may develop with dx2_y2 
order parameter at that special filling. The detailed study of the electron system at the Van Hove singularity shows 
that there is competition with magnetic instabilities, although for intermediate values of t' the pairing condensate 
develops at a higher energy scale and superconductivity prevails For the t — t' Hubbard model with t' < 0, we 
are led to propose a phase diagram with a superconducting instability with d^^-y^ order parameter above the Van 
Hove filling, for intermediate values of t' ( above ~ —0.3) . Below the Van Hove singularity we have found that there 
is a range of dopings where the superconducting instability has extended s-wave pairing symmetry, with nodes at the 
location of the inflection points. 

Our description of the anisotropic Fermi lines relies on the assumption that these do not change much upon lowering 
the cutoff Aq. The RG approach provides a consistent picture if only the electrons near the Fermi line are affected 
by the interaction, that is, Ao ^ J7 ^ £f, where ej;- is the Fermi energy. Different superconducting instabilities are 
known to arise in very dilute systems, for instance, where this approximation is not valid [S. 



V. SADDLE POINTS AT THE FERMI LINE 



In this section we analyze the instabilities of a two-dimensional electron system whose Fermi line passes through 
the saddle points (tt, 0) = A and (0, n) = B shown in Fig. 1 . To extract the physics of the saddle point we will take 
as an example the Hubbard t-t' model where perfect nesting is absent. We think that most of the features that we 
get are due to the presence of the saddle points in the Fermi line more than to the specific model. 

As mentioned in the previous section, this case can be seen as the limiting situation in which two inflection points 
merge into a saddle point; the technique developed previously can not be directly applied since the Fermi surface 
geometry becomes singular in this case. 

The presence of saddle points at the Fermi line has been related to the physics of the cuprates from the very 
beginning as a possible explanation for their high Tc superconductivity |3^j30| , |3^ , ^ , ^ , ^j2^ , [l^ , ^ , ^ . Its interest 
was reinforced by the photoemission experiments [^js) showing that the hole-doped materials tend to develop very 
flat bands near the Fermi level. The subject has evolved into the so-called Van Hove scenario that can be studied in 
the literature. In this section we will review the features of the model that can be extracted from a renormalization 
group point of view and derive its Kohn-Luttinger superconductivity following closely the scheme set in the previous 
sections. 

As described in section 2, the starting point of the RG study for fermion systems is the bare action despicted in 
(I): 

J dtd'p {i^+{p)dt'i'Ap) - {s{p) - eF)^+{p)^a{p)) 

+ J rftdVrf^P2rf^P3rfVC^(Pl,P2,P3,P4)*^(Pl)^'+.(P2)*a'(P4)*a(P3)'5(Pl +P2 -P3 -P4) , 

including the four fermion interaction term. 

Of the next two steps followed in section 2, which are the keypoints to obtain the Fermi liquid behavior, none can 
be implemented in our case. The first is associated to the isotropy of the Fermi line and consists in the kinematical 
decomposition of the momenta into parallel and orthogonal to the Fermi line at each point. It makes the scaling 
behavior of the meassure in the action effectively one-dimensional and is at the origin of the kinematical classification 
of the marginal channels. The second one is the Taylor expansion of the dispersion relation e(p) to the linear order 
what makes the momenta to have the same scaling behavior as the energy. In the present case the dispersion relation 
associated to the t-t' Hubbard model 
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e(k) = —2t [cos(fca:a) + cos(fcj,a)] — At' cos{kxa) cos(fcya) , (21) 

can be expanded around each of the saddle points A, B of Fig. 1 as 

sa.b(M) ^T{tT2t')kla^ ±{t±2t')kla^ , (22) 

where the momenta k^, ky measure small deviations from A, B. The quadratic dependence of the dispersion relation 
on the momenta induces the scaling tranformation: 

u —* suj , k — + s"'"/^k . 

Now it is easy to see that the same scaling law of the Fermi fields obtained in the Fermi liquid analysis: ^'(p) 
5~i/2\i>(p) , makes the free action marginal — since the scaling of the integration measure is the same — , and makes 
the four-fermion interaction to be marginal irrespective of the kinematics as now the delta function scales as the 
inverse of the momentum for a generic kinematics. 

Next, also in contradistinction with what happens in the marginal couplings of the Fermi liquid, the renormalization 
of the couplings in the Van Hove model is nontrivial due to the logarithmic divergence of the density of states dictated 
again by the dispersion relation. As described in section 3, the density of electron-hole pairs at a given momentum 
q depends on the curvature of the Fermi surface at the two points connected by q. 

This situation is similar to the one-dimensional model where there are also logarithmic singularities in the particle- 
hole and particle-particle propagators. As we will now see, the similarity stops at the classification level. The one loop 
analysis shows that singular quantum corrections arise in the saddle point model only when the momentum transfer 
in the process is of the forward (F), exchange (E) or BCS (V) type. Moreover the different channels do not mix at 
this level making the study more similar to the Fermi liquid case. 

In order to make clear the former statement we start by performing the standard analysis of the Van Hove g-ology. 
Although inclusion of spin is not necessary for the description of the Kohn-Luttinger superconductivity, we shall take 
it into account to include magnetic instabilities. 

The complete classification of the marginal interactions including the two flavors A and B follows exactly the one 
that occurs in the g-ology of one-dimensional systems [ ^0|j4l|| where the role of the two Fermi points is here played 
by the two singularities. 

In general there are four types of interactions that involve only low-energy modes. They are displayed in Fig. ^ 
where the interaction is represented by a wavy line to clarify the process it refers to. 

In the model that we are considering, with a constant interaction potential, the wavy lines should be shrunk to a 
point giving rise to couplings typical to the ^I*"* quantum field theory. The spin indices of the currents are opposite in 
all cases to stay as close as possible to the original Hubbard interaction. In this spirit, parallel couplings or spin- flip 
interactions shall not be allowed. 

All one loop corrections to a generic four fermion coupling are depicted in Fig. |[ 

There are four types of corrections in Fig. J^: direct (BCS) and exchange (EXCH) interactions (Fig. ^ (a), (b)), 
particle-hole interactions called RPA in Fig. ^(c), and vertex corrections called VERT in Fig. ^ (d). 

Once the interactions are shrunk to a point, they turn into the diagrams depicted in Fig. |lO| which are the ones to 
be computed in the one- loop calculation. We must keep in mind that, according to Fig. |^, the corrections induced 
by the BCS and RPA diagrams (Figs. ^ (a) and ^ (c)) have a minus sign relative to the others. 

All the one-loop diagrams can be parametrized in terms of the four polarizabilities involving particle-particle or 
particle-hole processes with the two lines corresponding to the same (zero momentum transfer) or to different flavors 
A, B (momentum transfer Q). These have been computed in |l^,^,Q, their cutoff-dependent part at energy a; = 
is 



2ttH °'e(q)' 
c' A 



2t:H ^'ta2(q-Q)2' 



c 



,2 



A 



Xpp(Q)-^arctan(2rf'A)log|^-5^A_| , (23) 



where 
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log 



(l + Vl-4(t'A)^) /{2t'lt) 



(24) 



Leaving aside for a moment the BCS graph, we see that all one-loop corrections are written in terms of the particle- 
hole polarizabilities (|2^) which diverge if the momentum transfer is zero or Q. That means that the RPA and VERT 
graphs will only provide corrections to a coupling if the kinematics of the vertex is such that pi ~ , P2 ~ P4 , 
i.e. both are forward processes in the Fermi liquid language which will renormalize only the amplitudes having the 
specified kinematics. The EXCH graph of Fig. ^will in turn be divergent only for the kinematics pi ^ Pi , P2 ^ Ps , 
i.e. it corresponds to an exchange process which renormalize exchange amplitudes. 

The analysis of the BCS channel is similar to the one done for the isotropic Fermi line. It is easy to see that 
particle-particle processes will provide a logarithmic renormalization only to those couplings whose kinematics is 
fixed to be of the BCS type. In the differential approach that we are using, this is best seen graphically as depicted 
in Fig. |ll|. 

Two energy integration slices contributing to the computation of the BCS graph are shown in Fig. |ll|. It is clearly 
seeen that, unless the total momentum of the incoming particles adds to zero, the area of the intercept of the two bands 
for which two high energy intermediate states contribute in the loop (which measures the cutoff dependence of the 
diagram) is of order (dA)^. This is different from what happens in one dimension, where the BCS graph contributes 
to the renormalization of all quartic couplings. Moreover, forward or exchange processes do not contribute to the 
BCS flow. 

We are then faced to a very similar situation as in the Fermi liquid case. Instead of talking of renormalization of 
couplings, we must analyze the fate of each channel forward, exchange and BCS, independently as they will not get 
mixed at the one loop level. 

Let us first analyze the F channel. It is easy to see that at the one loop level and for the case that we have chosen to 
analyze of point like interactions between currents of opposite spins, this coupling does not flow since we can not draw 
the diagrams of the corresponding corrections without invoking parallel couplings in the case of the RPA diagrams 
or spin flip interactions in the VERT diagram of ^. In the more general case in which parallel couplings are included 
from the beginning, there is an exact cancelation between the two types of F couplings as can be seen from Fig. |^. 
For each diagram of RPA type that can be drawn there exists another of VERT type that cancels it. This cancellation 
that, to our knowledge, was first noticed in the original work of ref. will occur to all orders in perturbation theory 
if the given conditions hold. Any k-dependence of the interaction would destroy this symmetry. In particular the 
very process of renormalization necessarily induces that dependence. In the presence of parallel couplings, the two 
graphs (c) and (d) in Fig. || will renormalize differently giving rise to a flow for the F channel. We will not address 
this question here. A related problem is discussed in 

The flow of the exchange channel must be examined by for each of the couplings of Fig. ^ that will from now on be 
denoted by E meaning that the momenta of the external legs are fixed to the exchange kinematics. The RG equations 
are obtained by "opening up" the graph and inserting the polarizabilities in such a way that the resulting graph is of 
the type of Fig. || (b) and the vertices are made up of the tree-level interactions of Fig. ^. 

Let us first discuss the behavior of any coupling, say i^intcr- ^'intcr is renormalized by the diagrams shown in Fig. 

Adding up the one-loop correction to the bare coupling we find the vertex function at this order 



2t:H 



inter 



Kmk) log 



(25) 



Following the usual RG procedure, we define the dressed coupling constant at this level in such a way that the vertex 
fimction be cutoff independent. Defining / = log A, we get the RG equation 



9ginter(A) 

dl 



2t:H 



(-^inter 



-^umk) ' 



(26) 



were the polarizability involved is the interparticle polarizability and the sign of the beta function is negative as 
corresponds to the "antiscreening" diagram of Fig. || (b). The same equation (^6|) is obtained by integrating over a 
differential energy slice as we mentioned earlier. 

By this method we obtain the following set of coupled differential equations for the E couplings: 



dE. 



1 



dl 

5i?back 



dl 



47r2 
1 

2^ 



-S'back) 



C (i^intra^^back) 
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Similiter ^ r' (f'^ 4- 



dl At:H 

c' (-Eintcr^^umk) (27) 



5i?umk 1 



where c, c' are the prefactors of the polarizabihties at zero and Q momentuni transfer, respectively given in (p^. 

Due to the relative plus sign of the corrections, all E couplings will grow if repulsive giving rise to instabilities in 
the system that can be identified by means of the response functions. It is easy to see that the E kinematics couples 
to the response functions corresponding to the order parameters 

^SDIV = X! Ks'^stO'k+Q,t , 
k 



The growth of these response functions drives the system to antifcrromagnetic and ferromagnetic ground states 
respectively. The phase diagram can be seen in ||2^ . 

Let us now see how the Kohn-Luttinger mechanism works in this case. As discussed in Section 2, superconducting 
instabilities will arise in the system by means of the BCS channel whenever a negative coupling develops. RG equations 
can be obtained for the BCS coupling by the same procedure followed with the E couplings. We will denote the BCS 
couplings by V implying the prescribed kinematics in the couplings of Fig. ^ and correct them at one loop by the 
BCS graphs in Fig. || As we are looking for pairing instabilities we will only consider in Fig. ^ (a) couplings with 
total momentum equals zero (not Q). This leaves us with Vintra and Vumk- Both involve the polarizability Xvvi^) 
( p3| ) . We will take care of the log^ dependence of the particle-particle susceptibility by the method discussed in . 
It consists in noting that the log^ can be factorized as the usual RG log and another one which is due to the divergent 
density of states exactly at the Fermi line and will take a constant value when the renormalization of the chemical 
potential is taken into account. We will denote this constant value by = log(A//i). Proceeding as before, we get 
for the V couplings the following set of equations: 

^Kmk 1 



dl 2ttH 



CK (MntraKmk) 



As there is no mixture of particle-hole channels, integrating this equation is equivalent to summing the ladder series. 

The flow of the couplings is shown in Fig. ^ It is clear that if all the couplings are set to a common value at the 
beginning, say the value U of the original Hubbard model, then l^ntia — Kimk all along and they would flow to zero 
along the diagonal in Fig. 

The Kohn-Luttinger mechanism will be at work if at some stage of the flow we get as initial condition Vintra ~ Vumk < 
0. This difference can be caused by the finite corrections induced in the BCS vertices by the particle-hole diagrams 
in Fig. ^ at the early stage of the flow. As discussed in that an operator is irrelevant means that it will flow to 
zero in the course of the renormalization, but it does not mean that it can be set to zero at the beginning without 
altering the physics. In this case, the BCS couplings l^ntra and Vumk will get finite particle-hole corrections from the 
exchange diagram that make them differ at intermediate values of the cutoff providing the initial conditions necessary 
for Kohn-Luttinger superconductivity. 

From equations (|2^) of the exchange couplings it is clear that under the initial condition that all couplings are 
equal to U, -Eintra = -Eback and -Eintcr = -^umk, all along the flow. As the corrections coming from the E vertex have 
a positive sign, we see that in the region of parameter space where c < c' in (p^, we will soon get Vintra < Vumk and 
then the RG flow of the BCS diagram will start with the initial condition necessary for d-wave superconductivity to 
set in. 

VI. CONCLUSIONS 

In this paper we have analyzed the instabilities induced on a strongly correlated electron system in two dimensions 
by the anisotropics of its Fermi line. Although we have chosen to exemplify the features with a Hubbard t — t' model 
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we believe that the study done is very general and will apply for other systems as well. Our main conclusion is that 
changes in the topology of the Fermi line of an electron system driven by doping, pressure or whatever means, give 
rise to instabilities of the system that will drive it away from Fermi liquid behavior. In the absence of nesting or 
other accidental symmetries, the instabilities will be predominantly superconducting. Although this idea is not new, 
as it goes back to the Kohn-Luttinger mechanism of the sixties, we have been able to see it at work in a microscopic 
model with different kinds of anisotropics. 

We have used a renormalization group approach which tracks the origin of the Kohn-Luttinger mechanism on the 
effect that some irrelevant operators have on marginal couplings. The scaling analysis allows us to predict, on very 
general grounds, a critical dimension of 3/2 for the onset of non-Fermi liquid behavior, a result that remains to be 
tested. Inflection points, a generic feature of some Fermi lines, give rise to a superconducting order parameter with 
extended s-wave symmetry with nodes located at the inflection points. 

The case of extreme anisotropy of the Fermi line having saddle points has been analyzed with the same technique, 
obtaining a very rich phase diagram for the t — t' Hubbard model with dx^-y^ superconductivity as the main instability 
in a certain range of parameters. The same symmetry of the order parameter is found for fillings above the saddle 
points, what suggests the possibility of a continuum transition in the overdoped regime of hole-doped cuprates. 

Despite the theoretical inspiration of this paper, the analysis performed contains some phenomenological predic- 
tions that can be tested experimentally in compounds which are well described by the t — t' Hubbard model. In 
particular, the closeness of a transition towards a d-wave and an extended s-wave superconductor implies that, in 
the absence of perfect tetragonal symmetry, a mixture of the two is likely. This possibility has been already discussed 
on phenomenological grounds |44[ . Different experiments, like photoemission |45| or interlayer tunneling in BSCCO 
p6| , can be interpreted in terms of the coexistence of different order parameters, or the existence of extended s-wave 
superconductivity. Note, finally, that a recent proposal for the Fermi surface of BSCCO, based on photoemission 
experiments |^], is close to the case with inflection points discussed here. Such a Fermi surface leads naturally to an 
extended s-wave order parameter in the superconducting phase. It would be interesting to see if this prediction is 
conflrmed by experiments. 



Note added — After this review has been published, the role played by the renormalized forward scattering 
interactions near a Van Hove singularity has been clarified in Ref. Furthermore, two detailed RG studies have 
appeared ^,|5^ supporting the idea that d-wave superconductivity takes place in the t — t' repulsive Hubbard model 
near the Van Hove filling. 
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CO 




FIG. 4. Available phase space for particle hole excitations with momentum Q = 2kF and energy less or equal to w. One of 
the two branches of the electronic band has been displaced by Q. 
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FIG. 5. Circular Fermi surface (top) and the two possible behaviors of Imx(Q,a;) at low energies. Middle: Generic case of 
non parallel patches of the Fermi surface. Shaded area is the available phase space. Bottom: Wavevector Q2 which connects 
two regions where the Fermi surface is parallel in the two patches. 
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FIG. 6. Top: Anisotropic Fermi surface with inflection points. Bottom: Available phase space when Q connects two 
inflection points. 
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FIG. 7. Diagrams building up the effective interaction in the BCS channel. 



19 




20 




A3 



RPA 

(c) 



Ai 




VERT 
(d) 




FIG. 9. Diagrams contributing to the one-loop order correction to tlie interaction potential. 



21 



A2 



A4 




A1 A3 
BCS 

(a) 



A4 A3 

>o< 

A1 A2 
EXCH 

(b) 



A3 



A4 




A1 A2 
RPA 
(c) 



A3 A4 

A1 A2 

VERT 
(d) 



FIG. 10. Diagrams of Fig. ^with the interaction shrunk to a point. 
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FIG. 11. Two energy slices used in the computation of the BCS graph. 
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FIG. 12. Renormalization of the different couplings described in the text. 
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FIG. 13. The flow of the V couplings described in the text. 
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